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Evidence for defect-mediated tunneling in hexagonal boron nitride-based junctions 
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We investigate tunneling in metal-insulator-metal junctions employing few atomic layers of hexag¬ 
onal boron nitride (hBN) as the insulating barrier. While the low-bias tunnel resistance increases 
nearly exponentially with barrier thickness, subtle features are seen in the cur rent-volt age curves, 
indicating marked influence of the intrinsic defects present in the hBN insulator on the tunnel¬ 
ing transport. In particular, single electron charging events are observed, which are more evident 
in thicker-barrier devices where direct tunneling is substantially low. Furthermore, we find that 
annealing the devices modifies the defect states and hence the tunneling signatures. 


Van der Waals heterostructures, where layered stacks 
of two dimensional materials are embedded in precisely 
desired patterns, have gained immense attention in re¬ 
cent times 1,2 . Such tailor-made structures of graphene, 
hexagonal boron nitride (hBN), metal dichalcogenides 
and other 2D materials offer exotic device geometries to 
explore new physics. Tunnel junctions with an atomically 
thin insulator barrier sandwiched between atomic layers 
of 2D materials form an interesting structure in this re¬ 
spect. In conventional two dimensional semiconductor 
double layer structures, tunneling has shown remark¬ 
able features, including resonant tunneling, Coulomb 
correlations at high magnetic fields and Landau-level 
spectroscopy 3 6 . In the regime of 2D layered materials, 
hBN with a band gap of ~ 5.9 eV 7 is an ideal candidate 
for an insulating barrier 8 . Recent studies on heterostruc¬ 
tures with single and bilayer graphene as electrodes and 
hBN as the insulator have shown interesting features, 
including a very strong negative differential resistance 
(NDR) 9-12 . In all of these transport studies hBN is as¬ 
sumed to be a benign element. However, from a materials 
perspective there have been extensive efforts to under¬ 
stand the underlying structure and defect mechanisms 
in thin layers of hBN, which can have important conse¬ 
quences on electrical transport 13-15 . Here we present de¬ 
tailed tunneling transport measurements on simple junc¬ 
tions consisting of metal or graphite electrodes separated 
by a thin hBN layer. Our results demonstrate that the 
hBN insulator can yield unexpected transport signatures 
suggestive of defect-mediated tunneling processes. 

In a conventional metal-insulator-metal (M-I-M) junc¬ 
tion, the tunnel current-voltage characteristic ( I-V ) is 
linear at low biases, with the tunneling resistance in¬ 
versely proportional to sample area and exponentially 
dependent on the barrier thickness d. Our studies show 
that the I-V with hBN as the barrier is distinct from 
such a simple behavior in various respects. Thinner bar¬ 
rier devices show a finite linear tunnel current at low 
biases and a roughly exponential dependence of the low 
bias resistance with d, complying with standard quan¬ 
tum mechanical tunneling, in agreement with previous 
reports 9,12 . However, in relatively thicker barrier devices, 



FIG. 1. (a) Device structure for h-BN tunnel devices with 
Cr/Au electrodes, (b) An optical micrograph of a thin hBN 
flake. Dotted line shows the thin region of interest, (c) hBN 
flake transferred on to a set of bottom electrodes (d) Tun¬ 
nel device with the horizontal top electrode forming multiple 
tunnel junctions, one of which is marked as a box. Scale bar 
is 10 /im. 


we find signatures of Coulomb blockade and single elec¬ 
tron tunneling events. In contrast to the conventional 
systems that show such effects, namely metallic islands 
and semiconductor quantum dots at low temperatures 16 , 
strikingly similar features are seen here in a presumably 
large area planar M-I-M junction. 

We have primarily investigated tunnel junctions com¬ 
prising Cr/Au-hBN-Cr/Au (see Fig. la) with 2-6 atomic 
layers of hBN forming the insulating layer, beyond which 
tunnel current was unobservably low. The device fab¬ 
rication begins by mechanically exfoliating hBN flakes 
on a Si/Si02 wafer. Figure lb shows a typical hBN 
flake. These are further characterized by a combina¬ 
tion of optical microscopy, atomic force microscopy and 
Raman spectroscopy to determine the number of lay- 
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ers. We note that employing a Nomarski interference 
contrast microscope helps in identifying thinner layers of 
hBN, which have lower optical contrast in comparison to 
graphene. The hBN flake is then transferred to another 
Si/Si02 wafer previously patterned with narrow Cr/Au 
electrodes, shown as vertical gold lines in Fig. lc. We 
have used the dry transfer method developed by Wang 
et. al 2 for all our devices. The top Cr/Au electrode 
(horizontal gold pad in Fig. Id) is then lithographically 
patterned. In all the devices, there are several bottom 
electrodes and one top electrode, thus forming multiple 
tunnel junctions (one out of four tunnel junctions in Fig. 
Id is indicated by the box) each of which can be measured 
in a four-terminal geometry. Typical junctions vary from 
1.4 x 0.5 fim 2 to 4 x 2 /im 2 in area. The tunneling cur¬ 
rent I and differential conductance dl/dV in response to 
voltage excitation (V = Vdc + V ac , with V ac = 0.5 mV 
at 13 Hz) were recorded simultaneously. Three differ¬ 
ent batches of hBN crystals were used for the present 
study and no qualitative differences were seen. In ad¬ 
dition to junctions with Cr/Au electrodes, we have also 
examined graphite-hBN-graphite devices. All of the data 
were taken at T — 4.2 K. 

We first focus on thinner hBN devices with <5 atomic 
layers (each layer being ^ 0.34 nm thick 8 .) Ten junc¬ 
tions (from four separate hBN flakes) with varying d and 
area A were tested. The I-V curves show an ohmic de¬ 
pendence at dc bias Vdc < ± 100 mV in all the junc¬ 
tions tested. (At higher biases, the I-V characteristics 
show substantial non-linearity. We will return to this 
later.) Figure 2a presents a optical-micrograph of a three 
atomic layer hBN device with junctions of varying area, 
along with I-V data from two of the junctions. At low 
bias voltages both junctions display a linear I-V. The 
lower inset to the figure shows that the current density 
J = I/A for the two junctions is the same in the linear 
regime \Vd c \ < 100 mV, as expected for a simple tunnel 
junction. 

At low biases with negligible barrier deformation, the 
linear current density for a M-I-M barrier 17 is, 

(-2d^/2rn*4^\ 

J oc exp ^--- )V, (1) 

where <pB is the tunnel barrier height and ra* is the effec¬ 
tive mass. Figure 2b reveals that the low-bias tunneling 
conductance per unit area G = J/V in our thinner bar¬ 
rier devices exhibits this basic exponential dependence. 
Not surprisingly, small variations in G among the various 
junctions with same d are observed, shown by error bars 
in Fig.2b. These variations might be due to undetected 
variations in the thickness d of the hBN layer or lack of 
uniform contact with the gold electrodes. An estimate 
to the factor m*(j)B can be obtained by fitting a simple 
exponential to the G vs. d data. Assuming the effec¬ 
tive mass to be 0.26mo for the conduction band edge in 
hBN 18 , we get <pB =3.3 eV. This seems reasonable given 
the presently unknown alignment of the ^5.9 eV hBN 
band gap relative to the Fermi level in the electrodes and 
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FIG. 2. (a) Tunnel current (/) as a function of dc bias (Vdc) 
for a typical thin hBN (dotted line) device with « 3 layers. 
Two tunnel junctions are shown in the inset optical micro¬ 
graph, with areas 3.8 x 2 ym 2 (red) and 2.3 x 2 ym 2 (blue). 
Scale bar is 10 ym. The inset graph shows the tunnel cur¬ 
rent density (J) for the two junctions as a function of bias, 
(b) Conductance per unit area G as a function of number of 
layers. Dotted line indicates an exponential fit. 


the appropriate effective mass for an electron tunneling 
far from the hBN conduction and valence band edges. 

We now turn to tunneling transport in samples with 
thicker hBN barrier layers (> 6 layers). Figure 3a shows 
a typical device with five tunnel junctions of area 1.4 x 
0.5 /rm 2 . Here, we find large discrepancy in the quantita¬ 
tive values of the tunnel current density between the dif¬ 
ferent junctions that are geometrically identical. While 
most of the junctions had a clear suppressed tunneling 
region around zero bias (green), some of the junctions 
also showed discrete steps in the I-V (red and blue). 
Figure 3b shows dl/dV as a function of V^ c . Consis¬ 
tent with the steps in I-V , periodic oscillations are seen 
with an energy spacing « 60 meV. This zero-bias sup¬ 
pression and staircase pattern in I-V bear strong resem¬ 
blance to Coulomb blockade features seen in metallic or 
semiconducting dots 16 . While the suppression of current 
around zero bias is a universal feature of Coulomb block¬ 
ade, the Coulomb staircase constitutes a rather special 
case where the tunnel barriers that isolate the island from 
the leads are asymmetric (with unequal tunneling resis- 
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FIG. 3. (a) Tunnel current (/) and (b) differential conduc¬ 
tance (dl/dV) as a function of dc bias (Vdc) for a typical thick 
hBN device with « 6 layers. Green, blue and red denote three 
junctions with area 1.4 x 0.5 fim 2 as shown in the inset. Scale 
bar is 5 fim. 


tances and/or capacitances) 19,20 . In such a scenario, an 
extra electron dwells in the island as the transmittance 
through one of the barriers is suppressed. Coulomb stair¬ 
case patterns are observed with voltage spacing given by 
e/C, where C = C\ + C 2 + C s , C\ and C 2 being the 
capacitances to the source and drain leads and C s is the 
self capacitance of the island. The tunneling character¬ 
istics showed very little dependence on temperature or 
magnetic field, with the Coulomb blockade features dis¬ 
appearing above about 77 K. 

Sharp features in dl/dV arising from various inelas¬ 
tic processes involving phonons have long been used as 
a sensitive tool to probe the phonon spectrum of the in¬ 
sulator, the electrodes or molecules trapped at the inter¬ 
face 21 . However, the periodicity that we find in the struc¬ 
tures, along with the strong blockaded region around 
zero bias suggest that these are most likely signatures 
of single electron charging events. The charging energy 
E c = e 2 /2 C is typically between 20 — 100 meV in most 
samples, giving a capacitance ~ 0.8 — 4 aF, which indi¬ 
cates a very small region as the main contributor to the 


transport, even though the tunnel junction itself is of the 
order of microns. Ignoring the capacitance to the leads 
for now, the self capacitance of a disc of radius r is given 
by C s = 8eoe r r, which for a capacitance of 1 aF and 
e r = 4, would give a radius of 3.5 nm. The question to 
ask then, is what component here acts as a dot or an is¬ 
land separated from the metal electrodes? Interestingly, 
one possible candidate that complies with the small size 
of the island are the intrinsic defects in the hBN layer. 
hBN crystals have been known to contain carbon and 
oxygen impurities in the order of 10 17 /cm 3 , which have 
been experimentally determined by secondary ion mass 
spectroscopy and cathode luminescence 22 . This would 
correspond to about 100 such defects in a lxl /im 2 area 
for a 1 nm thin sample of hBN. In addition, transmis¬ 
sion electron microscopy (TEM) studies have revealed 
vacancies, interstitial defects and ionized centers in hBN 
samples 14,23,24 . Recent reports of photo-induced doping 
in graphene on hBN substrates 25 and scanning tunneling 
microscopy studies of defects in exfoliated hBN 26 sub¬ 
stantiates this picture. Importantly, defects and voids in 
hBN can even lead to interlayer bonding in contrast to 
few layer graphene sheets 15 . These suggest that defects 
can play a major role in tunneling transport as well, form¬ 
ing small isolated islands in the insulating matrix acting 
like quantum dots, where adding an extra electron costs 
energy. We thus have two scenarios: (1) conventional 
direct tunneling that scales with area and depends ex¬ 
ponentially on the barrier thickness and (2) tunneling 
mediated by defects. The former is more pronounced 
in thinner samples, while the latter contributes to the 
sharp features in dl/dV in thicker samples where direct 
tunneling is suppressed. Note that the capacitance to 
the electrodes is a relevant quantity here, as the hBN 
barrier is very thin with a thickness of ~ 1.8 — 2 nm. 
The capacitance to the leads for a disk of radius 3 nm 
is given by C\ = C 2 ~ Aeoe r /d = 1 aF, for d = 1 nm 
while, the self capacitance C s = 0.8 aF. Indeed, the metal 
electrodes in our layered tunnel junctions very effectively 
screen the hBN from externally applied electric fields. 
This has the unfortunate consequence that the “Coulomb 
diamonds” which are characteristic of single electron tun¬ 
neling events, cannot be observed in our devices. 

In order to eliminate any possible influence of the 
electrode metal or the fabrication process, we tested a 
graphite-hBN-graphite device consisting of three junc¬ 
tions with « 6 atomic layers of hBN. Bottom graphite 
electrodes were fabricated by mechanical exfoliation fol¬ 
lowed by e-beam lithography and oxygen plasma etch to 
create multiple electrodes. The hBN flake was deposited 
on the bottom graphite electrodes, followed by the trans¬ 
fer and fabrication of a top graphite electrode. Graphite 
offers an atomically flat electrode surface for the tunnel 
junction. As can be seen from Fig. 4a, we find very sim¬ 
ilar characteristics even in this case, with a blockaded 
region around zero bias and step-like characteristics in 
the I-V. Thus, a bad contact with the Au electrodes or 
the possibility of a Au dot in the hBN layers due to the 
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FIG. 4. (a) Tunneling I-V for a typical graphite-hBN- 

graphite junction with an area 2 x 1.5 fi m 2 . Two other junc¬ 
tions showed similar characteristics, (b) I-V for a thick hBN 
barrier tunnel junction with Cr/Au electrodes, before (red) 
and after (blue) annealing in Ar-hG environment. Similar re¬ 
sults were obtained for three more junctions in two sets of 
devices. 


fabrication procedure can be safely dismissed. 

In addition to probing the signatures of defects in 
charge transport, we also find that we can manipulate 
them by external means. To achieve this, we annealed 
some of the devices that showed significant Coulomb 
blockade features in an Ar:H 2 environment at 350° C for 
3-4 hours. Figure 4b shows the results for one such junc¬ 
tion; similar results were obtained from other junctions 
as well. As the figure reveals, the pristine sample, for 
which the hBN layer was sufficiently thick that no di¬ 
rect tunneling was expected or observed, showed pro¬ 
nounced Coulomb blockade and staircase features, with 
E c ~ 100 meV. After annealing all signatures of single 
electron charging effects are absent. Only at relatively 
large bias, \Vd c \ ~ 0.6 V, does a significant tunneling 
current gradually appear. In graphene and other lay¬ 
ered 2D materials, annealing in an atmosphere of Ar:H 2 
or Ar :02 is expected to remove leftover polymers and 
other residues from device fabrication and improve ad¬ 
hesion to the substrate 27 . Annealing the exfoliated hBN 
before transferring graphene on it is often considered an 
important step in achieving high mobilities for graphene 



FIG. 5. Differential conductance dl/dV vs. Vdc for the thin 
barrier device shown in Fig. 2a. Red and blue denote the two 
junctions shown in the inset to Fig. 2a. 


devices 28 . However, the microscopic mechanisms behind 
such drastic improvements with annealing are still not 
fully clear. In our devices, annealing seems to elimi¬ 
nate the defects in the active tunnel junction area. In 
bulk hBN, vacancy migration is expected to occur around 
500°C 29 . Recent TEM studies on atomic layers of hBN 
have revealed grain boundaries, square-octagon (4|8) de¬ 
fects, pentagon-heptagon defects (517) in addition to va¬ 
cancies and interstitials 13,14 . While experiments have re¬ 
vealed that 4|8 defects are mobile at about 800°C under 
e-beam irradiation 14 , 5|7 defects are expected to be move 
to the boundary and stop or vanish 30 . Whether anneal¬ 
ing cures some of the defects or mobilizes them will need 
a microscopic analysis using techniques like TEM and is 
beyond the scope of this work. Nevertheless, this indi¬ 
cates that the features have an intrinsic origin and can 
also be tuned via simple annealing steps. In addition, 
these are the first signatures of a direct influence of hBN 
on electrical transport in Van der Waals heterostructures. 

We now return to the thin barrier devices at large bias 
voltages. Beyond the ohmic low bias region, the I-V 
curves showed a pronounced non-monotonic region, in 
many cases resembling an NDR feature as seen in Fig. 
2a. Figure 5 shows the differential conductance dl/dV 
as a function of Vd c for the two junctions discussed in 
Fig. 2a, where the non-ohmic behavior can be seen be¬ 
yond Vdc ~ 0.2 — 0.3 V. While the red curve shows a 
clear NDR region in the negative voltage side, the effect 
is less pronounced in the positive side and for the blue 
curve. The voltage range probed is believed to be below 
the Fowler-Nordheim regime, where the barrier is essen¬ 
tially triangular, due to a very high applied bias 8 . In 
recent studies, NDR signatures were observed in single 
and bilayer graphene based heterostructures, which were 
attributed to resonant tunneling via momentum conser¬ 
vation when the energy bands of the top and bottom 
graphene layers were aligned 11,12 . However, we find sim- 
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ilar I-V curves in simple M-I-M junctions here, albeit 
with peak to valley ratios lower than in graphene based 
devices. Interestingly, the only common feature in all 
these devices is the hBN barrier layer. The origin of 
these non-ohmic features remains ambiguous. Finally, 
we also observe several weak peaks at regular intervals 
in the dl/dV (marked with arrows), in many of the sam¬ 
ples. These may be due to defect-mediated single electron 
charging events adding to the much larger direct tunnel 
current. 

In conclusion, we have reported tunneling transport 
experiments on simple tunnel junctions consisting of 
metal (or graphite) electrodes separated by hexagonal 
boron nitride barrier layers. For thin hBN layers, di¬ 
rect tunneling, with reasonable area and barrier width 
scaling is observed at low biases, while unexplained non¬ 
linear effects are seen at higher biases. In junctions with 


thicker hBN barriers, direct tunneling is negligible but 
strong evidence for defect-mediated single electron tun¬ 
neling processes is found. The defect states, and hence 
the tunneling transport in these thicker barrier devices, 
can be altered through a simple annealing process. 
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